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Abstract 

Given a planar curve defined by means of a real rational parametrization, we 
prove that the affine values of the parameter generating the real singularities 
of the offset are real roots of a univariate polynomial that can be derived 
from the parametrization of the original curve, without computing or making 
use of the implicit equation of the offset. By using this result, a finite set 
containing all the real singularities of the offset, and in particular all the 
real self-intersections of the offset, can be computed. We also report on 
experiments carried out in the computer algebra system Maple, showing the 
efficiency of the algorithm for moderate degrees. 


1. Introduction 

Intuitively, the offset curves to a given curve are “parallel” curves to the 
original curve, called the generator curve, at a certain distance. The offset¬ 
ting operation is important in Computer Aided Geometric Design (CAGD), 
because it can be used to give “thickness” to an object, in this case a curve, 
which is thin; the offsetting distance can be regarded as the desired “thick¬ 
ness” of the new object. Furthermore, offsets also have applications in fields 
like robotics or manufacturing [16j . |25j. 
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When we compute the offset of a curve, typically we want to reproduce in 
the offset the shape of the generator. However, sometimes the offsetting oper¬ 
ation introduces singularities that destroy the topology of the original curve. 
Hence, one has to identify the parts of the offset that should be discarded, 
and trim them away in a post-processing step. For example, in Figure [I] (left) 
we can see the parabola y = x 2 (in thin line) and its offset at distance d = 1 
(in thick line). We observe that while one of the connected components of the 
offset has the topology of the parabola, the other component has a different 
topology. This last component has one self-intersection and two cusps. Af¬ 
ter trimming away the loop containing the three singularities (see Figure [lj 
right), the two components of the new curve that we get, called the trimmed 
offset , have the topology of the original curve. Furthermore, in Section 4 of 
m it is proven that the trimmed offset can be easily computed whenever 
the parameter values of the self-intersections of the offset, or even a finite set 
containing them, are known. So the computation of the self-intersections of 
the offset is strongly related with the trimming operation. 

Trimming is important in computer aided design to keep the original 
shape and therefore improve the appearance of the image, but also in man¬ 
ufacturing. More precisely, if a certain curve is to be machined by the cylin¬ 
drical cutter of a milling machine, the cutter follows a certain-line trajectory 
specified by the offset, where the offsetting distance equals the cutter radius 
(see page 162 of [IS])- However, if the offset has self-intersections then the 
offset will have loops (like the one in the offset to the parabola, see Figure [lj 
left), giving rise to the problem of “gouging”, as it is called in NC machining 
(see Section 11 in [25]). Essentially, gouging implies that these small loops 
must be removed later by using a smaller size cutter. 

In this paper we deal with the problem of computing the singularities of 
the offset to a generator curve given by means of a rational parametrization, 
as it is common in CAGD. More precisely, we are interested in computing the 
affine values of the parameter giving rise to the offset singularities. In this 
context, a first difficulty is the fact that the offset does not need to be rational; 
in fact, if the offset is rational then the computation of its singularities is 
relatively easy, and can be done for instance by using the method in [26] . 
Since the offset to an algebraic curve is also algebraic, one might try to 
overcome the aforementioned difficulty by working with the implicit equation 
of the offset, in order to derive the singularities of the offset from there. 
Nevertheless, the offsetting process causes kind of an “explosion”, so that the 
implicit equation of the offset is much more complicated than the original 
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Figure 1: Offset (left), and trimmed offset (right) of the parabola at the 
distance d — 1. 


one, even for very simple curves. Therefore, deriving the singularities of 
the offset from its implicit equation is very time-consuming and quite often 
impossible in practice. 

As a consequence, many papers in the literature have proposed methods 
to derive the singularities of the offset by computing them directly from the 
generator curve, or by approximating the offset by another object which is 
easier to manipulate. In this sense, one should distinguish between “local” 
singularities, i.e. singularities, like for instance cusps, that are due to local 
phenomena, and self-intersections, which are due to the intersection of dif¬ 
ferent branches of the offset. Local singularities of offsets to rational curves 
are easy to find, since the values of the parameter t generating them are the 
solutions of the equation k(t ) = —1/d, where d is the offsetting distance and 
k(t) is the curvature of the curve. For singularities of this type coming from 
regular points of the generator, one can see Section 2.5 of mi. In fact, the 
results in p3] are applicable not only to rational curves, but to parametric 
curves defined by means of a regular parametrization, i.e. a parametriza- 
tion where the speed vector does not vanish. For the analysis of local offset 
singularities coming from singular points of the generator of an algebraic 
curve, one can see [2J. However, finding the self-intersections of the offset is 
a difficult problem. 

In [15] and [21], the self-intersections are directly derived from the gen¬ 
erator. In [15], this is done for the case where the initial curve admits a 
polynomial parametrization. In order to do this, the parameter values gener¬ 
ating the self-intersections are proven to be the roots of a polynomial which 
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is the quotient of a big determinant and a product of two polynomials. How¬ 
ever, the generalization to the case when the generator admits a rational, 
non-polynomial, parametrization is cumbersome. 

In 12T], the self-intersections of the offset are computed by solving a sys¬ 
tem in four variables and four equations, also derived from the generator 
curve. The idea is also applicable to a general regular parametrization, not 
necessarily rational. Nevertheless, in that case one has to deal with the 
numerical problem of approaching all the solutions of a nonlinear system. 
Notice that if a non-rational parametrization is used, the system is not nec¬ 
essarily algebraic. However, if the parametrization of the original curve is 
polynomial, the system is certainly algebraic. In this case, by writing the 
parametrization of the original form in Bezier form, in turn one can write 
the equations of the system as Bernstein polynomial equations. Then the 
solution of the system is reduced to finding the intersection of two bivariate 
Bezier patches with a certain plane. In order to do this, de Casteljau subdi¬ 
vision methods coupled with rounded interval arithmetic are involved. The 
generalization to rational curves, though stated to be feasible, is mentioned 
as a topic of future research. 

A second possibility, that has been explored by many authors, is to ap¬ 
proximate the offset by means of a simpler object, and then approximate the 
self-intersections, and therefore the trimmed offset itself, from that object. In 
some cases the offset is approximated by a rational curve [33], or a polynomi- 
ally parametrized curve m In other cases [10], (SO. 1231 an approximation 
with a polygonal line is used. In [19], the input is a planar rational curve, 
and a Gd-continuous biarc approximation of the curve is employed. Some 
other approaches to the problem and additional references can be found in 
Section 11.2.4 of [22] , 

In this paper, we provide a new method to find a finite set containing all 
the affine parameter values giving rise to real, non-isolated singularities of 
the offset. Our method computes these values from the generator curve, and 
does not require to compute or make use of the implicit equation of the offset. 
We were inspired by the ideas of HU, which in turn is related with [6], where 
the computation of the genus of the offset from the genus of the generator 
curve is addressed. The main idea of the method is the following: by na, 
one can establish a birational mapping between the offset and a much simpler 
curve. However, at the self-intersections of the offset, this mapping cannot 
be inverted. Additionally, one can prove that the same holds not only for 
the self-intersections of the offset, but for all real singularities of the offset. 
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From a computational point of view, our algorithm uses subresultants as an 
essential tool, jointly with root finding. Furthermore, we have implemented 
and tested our algorithm in the computer algebra system Maple; the code 
can be freely downloaded from [12] . 

Compared to other methods, our algorithm does not use any approxima¬ 
tion of the offset, and can be applied to possibly singular, rational curves. The 
algorithm has at least two advantages: first, the implementation is easy and 
requires only a few lines of code. Second, the description is, unlike ra. bd. 
basically the same regardless of whether the parametrization is polynomial 
or rational. In fact, the algorithm is presented under the assumption that the 
parametrization is non-polynomial. As a disadvantage, we can mention the 
potential appearance in the output of superfluous values of the parameter. 
This does not happen when the parametrization is polynomial, but it can 
happen, in certain cases, when it is non-polynomial. Nevertheless, even in 
this case our results can be applied to the offset trimming problem, since the 
presence of superfluous values does not affect the final result. We must also 
observe that when approached from a symbolic or symbolic-numeric point 
of view, the problem is inherently difficult, because the degree of the offset 
of a curve of degree n is bounded by 2(3 n — 2) (see Theorem 3.6 in [15] ). 
and the growing of the coefficients can be serious. For instance, the offset 
of the Descartes’ Folium x 3 + y 3 — 3 xy = 0 generically has degree 14 and 
114 terms. In spite of this fact, our experiments show a good performance of 
our algorithm for moderate, but far from trivial, examples. These examples 
include some curves analyzed in {15j, [21], in order to compare with those 
methods. 

The structure of the paper is the following. Generalities on offsets and 
subresultants are provided in Section [2] The strategy behind our method is 
presented in Section [3j This strategy is at first aimed to compute the self¬ 
intersections of the curve; however, we prove that, as a by-product, we also 
get the remaining singularities of the offset. The main result of the paper, 
jointly with the algorithm it gives rise to, are given in Section [4j details on 
examples, as well as a thorough analysis of the complexity of the algorithm 
and the growing of the coefficients, are also provided here. The conclusions 
of the paper are presented in Section [5] Although our algorithm can be easily 
described and implemented, the proof of the result it is based on takes certain 
work; the parts of the proof that are not essential to understand the main 
result of the paper and the subsequent algorithm, are given in Appendix I, 
Appendix II and Appendix III. 
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2. Preliminaries and generalities. 

2.1. The offset curve 

Let C be a real, rational plane curve, parametrized by 



( 1 ) 


where X(t),Y(t),W(t) are polynomials with real coefficients, and 


gcd(X(t),Y(t),W(t)) = l. 


We also assume that <f(t) is proper, i.e. birational or equivalently injective 
except perhaps for finitely many values of the parameter t. This condition 
ensures that C is a reduced curve (see Theorem 4.41 in (32]). Note that 
properness can always be achieved by reparametrizing the curve, if necessary 



The offset to C at distance d G M + , Od(C), is defined as the Zariski closure 
of the set of points (x,y) = 4>d(t), where 



Mt) = 


with 


U(t) = X'{t)W{t) - X(t)W'(t), V(t) = Y\t)W(t) - Y{t)W'(t). (3) 

Furthermore, in the paper we will assume that t is a real value. Therefore, 
(x, y ) = (j>d(t) means that the Euclidean distance between (x, y) G Od(C) and 
the point p = <f{t) G C, measured along the normal line to C through p = <f(t), 
is d\ we say then that p = <f(t) generates (x,y). When the first sign of ± 
and =F in the expression (|2]) is considered, the geometrical locus described is 
called the exterior offset ; if the second sign is chosen, the geometrical locus 
described this way is called the interior offset. Furthermore, 4>d(t) can be 
extended to the t 0 values where U 2 (t 0 ) + V 2 (t 0 ) = 0 by just taking limits 
t —^ t 0 (in the usual topology); the points computed this way also belong to 
the offset. Additionally, if lim t _ ) . oo 0(t) is an affine point, which happens iff 
deg(X(t)) < deg(W(t)) and deg(P(t)) < deg(W(t)), lim t ^oc4>d(t) generates 
two more points, which we will denote as P±oo, also belonging to the offset. 
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The computation of an implicit equation F(x, y) = 0 of Od(C) is addressed 
in |15j . Let us review some of the ideas of ca, which are relevant for our 
purposes. In [15] the following polynomials are introduced: 

P(x,y,t) := U(t)(W(t)x - X(t)) + V(t)(W(t)y - Y(t )) = 0, (4) 

Q(x, y , t) := (W(t)x - X(t)) 2 + (W(t)y - Y(t)) 2 - d 2 W 2 (t) = 0. (5) 

For the values of t which satisfy W{t) ^ 0, ( U(t ), V(t)) ^ (0, 0), the equation 
Q represents the normal line to C at the point Wpj) > w hil c the equa¬ 
tion (J5]) represents the circle of radius d centered at the point w§)\ 

The implicit equation of 0^(0) is determined by eliminating the variable t 
in the system formed by Q and (|5]). However, in order to avoid extraneous 
components (see ra for details), we must divide first P(x,y,t), Q(x,y,t) by 
their contents with respect to tjj Since Q{x,y,t ) can be written as 

Q(x,y,t ) = W 2 (t)(x 2 + y 2 ) — 2W{t)X{t)x — 2W{t)Y(t)y 
+X 2 (t) + Y 2 {t) — d 2 W 2 (t), 

one can easily see (Lemma 3.1, [15]) that the t-content of Q is equal to 

y(t)=gcd(W(t),X 2 (t)+Y 2 (t)). (6) 

Similarly (Lemma 3.2, [T5]) the f-content of P is equal to 

P(t) = (7) 


with 


a{t) = gcd{W{t),W'{t)), 7(t) = gcd{U{t)/a{t),V{t)/a{t)). (8) 


Let P(x,y,t ) and Q(x,y,t ) be the polynomials obtained after removing the 
f-contents from P and Q, 


P(x,y,t ) 


P(x,V,t) 

m 


Q{x,y,t) 


Q(x,y,t) 

y{t) 


(9) 


4 Let f(xi ,..., x r , x r + 1 ,..., x n ) be a polynomial in the variables xi,..., x r , x r +i ,..., x n 
with coefficients in a unique factorization domain. The content cont^...^ (/) of / with 
respect to xi,..., x r is the gcd of the coefficients of /, seen as a polynomial in x r +i,..., x n 
whose coefficients are polynomials in x\,...,x r . The polynomial / = ^ • / is 

called the primitive part of / with respect to X\, , x r . 
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and let 


H(x, y ) = Res t {P(x, y, t ), Q(x, y, t)). 

If pit) is constant, then F(x,y ) = H(x,y). If deg(//(i)) > 0 then H(x,y) can 
have extraneous, linear factors. Hence F(x, y ) is the result of dividing H(x, y ) 
by the product of the extraneous factors. Notice that if <f>(t) is polynomial , 
i.e. if W(i) = 1, then p{t) is constant and no extraneous factors can appear. 

In order to see that F(x, y) — 0 really corresponds to Od(C), one observes 
that F(x,y ) is a factor of Res t (P(x,y,t),Q(x,y,t)). By well-known proper¬ 
ties of resultants (see HU ) then for each (x 0 ,2/o) E Od(C) either there exists 
t 0 such that P(x 0 ,y 0 ,t 0 ) = Q(x 0 ,y 0 ,t 0 ) = 0, or (x 0 ,y 0 ) is a common zero of 
the leading coefficients lc t (P), 1 q(Q) of P(x,y,t ) and Q(x,y,t ) with respect 
to t. The next result, proved in Appendix I so as not to stop the flow of 
the paper, characterizes the points where both lc t (P), lct(<5) simultaneously 
vanish. 

Lemma 1. The only points (xo, yo) where the leading coefficients of P(x, y, t), 
Q(x,y,t ) with respect to t simultaneously vanish are P±oo, in the case when 
P± oo are affine points. 

Therefore for every (xo,yo) E Od(C), with perhaps the exception of P±oo, 
there exists t 0 such that P(x 0 ,y 0 ,t 0 ) = Q(x 0j y 0 ,t 0 ) = 0. In turn, this means 
that (x 0 ,2/o) simultaneously belongs to the normal line to C through p 0 = 
4>{t 0 ), and to the circle centered at p 0 of radius d, which implies that {x 0 ,y 0 ) 
is generated by t 0 via 4>d(t). Note that in that case W(t 0 ) ^ 0: indeed, since 
by hypothesis gcd(A"(f), Y(t), W{t)) = 1, W(t 0 ) = 0 implies that either X{t) 
or y(t) tend to infinity as t —» t 0 . 

A more algebraic definition of the offset curve, using an incidence dia¬ 
gram, can be found in [28]. We also summarize some notions and results of 
|5] and (28j, that we will use in this paper. First, Od(C) has at most two com¬ 
ponents. An irreducible component of 0 ( i(C ) is said to be simple if almost 
every point of that component is generated by just one point of C; otherwise, 
the component is called special. If C is properly and rationally parametrized, 
then Od(C) is reducible iff U 2 [t) + V 2 (t), in our notation, is a perfect square 
(see Corollary 3.4 in [5]). Furthermore, if Od(C) is reducible, then it has two 
rational components. Finally, if Od(C ) is irreducible then it is simple. 

If Od(C) is reducible the problem which we address in the paper can be 
solved in an easier way. Indeed, denoting the rational parametrizations of 
the components of Od(C) by aq(t) and x 2 (t), one just needs to: (1) com¬ 
pute the singularities of aq(f) and x 2 [t) separately (see for instance [29]); 
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(2) compute the intersections of the components of Od(C), for instance by- 
setting X\{t) = X 2 (s) and applying elimination methods. Thus, in the rest 
of this paper we will assume that Gd(C) is irreducible, in which case it is also 
simple. Therefore, the real part of Od(C), which is the object that we want 
to study here, corresponds to the points generated by f G R (see Remark [I] 
in Subsection 2.3), possibly with the exception of some isolated singularities 
and P ±oc . 

In this paper we will address the computation of the real affine, non¬ 
isolated singularities of Od(C), which are generated by real values of the 
parameter of 0(f) via ([2]). The only exception to this are the points P± OOJ 
which may be generated only by f = oo. Since our goal is to find the f values 
generating the singularities of Od(C), this is not really an issue. Nevertheless, 
if one also wants to detect whether or not P ±^ are singular, it suffices to 
reparametrize the curve so that P ±are generated by affine values of f, and 
then examine these f values. 


2.2. Subresultants. 

We refer to m, 0 , usi and [35] for further reading on the notions and 
results in this subsection. Let D be an integral domain, and let /, g G D [f] 
be the polynomials 

/(f) = a n t n + a n _if" 1 + • • • + ao, g(t) = b m t ,n + & m _if m 1 + • • • + bo, 


where deg(/) < n, deg (g) < m. 

Definition 2. For i G {0,..., inf (n, m) — 1} ; the Sylvester matrix of index 
i associated to /(f), n, g(t) and m, denoted by SylVj(P, n, Q, m), is the 
( n + m — 2 i) x (n + m — i) matrix: 


Syl v i(f,n,g,m) 



( a n 

.. ao 

\ 

1 


d n • • 

• «0 


bm 

.. b 0 


1 

\ 

bm • • 

• bo) 

I 


The Sylvester matrix of index 0 associated to /(f), n, g(t ) and m is denoted 
by Sylv(/, n, g, m). 
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If deg(/) = n and deg(g) = m then the Sylvester matrix of index 0 is 
simply called the Sylvester matrix of fit) and g(t), denoted by Sylv(/, g), 
and the Sylvester matrix of index i 0 is denoted by Sylvff, g). 

Definition 3. The determinant o/Sylv(/, g), denoted by Res (f,g), is known 
as the resultant of f(t) and g{t). 

The concept of determinant polynomial associated to a matrix provides 
one of the usual ways to define subresultant polynomials. 

Definition 4. Let A be a m x n matrix with m < n. The determinant 
■polynomial of A, detpol(A), is defined as: 

n—m 

detpol(A) = det(A k )t n ~ m ~ k 

k =0 

where A^ is the square submatrix of A consisting of the first m — 1 columns 
and the (k + m)-th column. 

Definition 5. The i-th subresultant polynomial of fit) and gift), Subres,(/, n, g, m), 
is the determinant polynomial of SylVj(/, n, g, m), i.e. 

SubreSj(/, n, g, m) = detpol(SylVj(/, n, g, m)), for 0 < i < inf (n,m) — 1. 

The sequence {Subresj(/, n, g, m)}*>o is called the subresultant chain of f,g. 

We have that deg(Subres i(f,n,g,m)) < i. The coefficient of degree i of 
the polynomial Subres ff, n, g, m), denoted by sres*(/, n, g, m), is called the 
i-th principal subresultant coefficient of /, n and g, m. When sres *(/, n, g, m) = 

0, the polynomial Subres i(f,n,g,m) is said to be defective. 

Theorem 6. [Fundamental Property of subresultants] Let K be the fraction 
field of D, and assume that deg(jf) = n or deg(p) = m, gcd (f,g) ^ f and 
gcd (f,g) ^ g. Then the first subresultant polynomial different from zero 
in the sequence {Subres k(f,n,g,m)}k>o is non-defective, and equal to the 
greatest common divisor of f,g in K[t]. 

Theorem 7. Let if : A —$■ B be a ring homomorphism and P,Q e A[t] be 
two polynomials with deg(P) = p and deg(Q) = q. If deg(ip(P)) = p and 
deg (p[(Q)) = q then for any i < q we have Subres i ('ip(P),p,'ip(Q),q) = 
■0(SubreSj(P,p, Q, q)). If deg(ip(P)) = p arid deg (ip(Q)) = q* < q (or 
vice versa), Subres i(t[(P),p,if(Q),q*) and ^(Subres^PjP, Q, q)) are pro¬ 
portional. 
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2.3. Singularities. 

Let V be an algebraic planar curve, implicitly defined by f(x,y) = 0. 
An affine point S G V is a singularity iff f x (S ) = f y (S ) = 0; a nonsingular 
point of V is said to be a regular point. One can always compute a local 
parametrization (x(h), y(h)) of T> around any point S G T>, regular or singu¬ 
lar, where x(h),y(h) are two analytic functions in a neighborhood of h = 0. 
An equivalence class of irreducible local parametrizations around S is called a 
place pfB]; we say that the place is centered at S = (x(0), y(0)). Furthermore, 
we say that the place is real if there is a representative of the class where all 
the coefficients are real. If (x(h),y(h)), with x(h),y(h) real, represents a real 
place, where x(h),y(h ) converge for \h\ < e, the set of points of D defined by 
(. x(h),y(h )) for \h\ < e, is called a real branch of V through (x(0), y(0)). 

If S' is a self-intersection of T> , then S is the center of several, different, 
places of T>. If I? is not parallel to the y- axis, one can prove [2] that any place 
centered at S can be written, in a coordinate system centered at S where the 
x-axis coincides with the tangent to V at S (see Figure |2j left), as 

V{h) = {h?,p q h* + ••■), 

with p, q e N, p > 1, q > p. We will say that V(h) is singular if p > 2. A 
point S £ TI is a singularity iff it is either the center of one singular place, 
or the center of several singular or regular places. In the first case, we will 
say that S is a local singularity of T>. In the second case, we will say that 
S' is a self-intersection of T>. Notice that S can simultaneously be a local 
singularity and a self-intersection when it is the center of several places, and 
at least one of them is singular (see Figure [2j bottom-right). 

Remark 1. In the rest of the paper we will address the computation of the 
real local singularities and the real self-intersections that are not isolated, i. e. 
that are the center of at least one real place. Since we are assuming that 
Od{C) is simple, then any real branch of Od{C ) comes from a real branch of 
C /If. Since in turn every real branch of C is generated by real values oft via 
4>(t), any affine t value generating a real, non-isolated singularity of Od(C) 
is real. 

3. Strategy for computing the offset singularities. 

3.1. The idea. 

Let us consider first the computation of the real self-intersections of 
OdiC). Later on we will see that our method to compute these singulari- 
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Figure 2: Local singularities and self-intersections. 


ties provides also, as a by-product, the local, real singularities of Od(C). Our 
idea to find the self-intersections of Od(C) is inspired in the strategy used in 
[17] to find the genus of Od(C). In [17], the curve implicitly defined in the 
(t, a) plane by 

a 2 - (U 2 (t) + V 2 (t)) = 0, 

is introduced. In order to avoid difficulties in the analysis of the Lvalues 
generating local singularities of C, we will use instead the curve At defined 
as 

a 2 — (U 2 (t) + V 2 (t)) = 0, 


with 


U = 


U 


and V = 


V 


gcd (U,V) gcd (U,vy 

Notice here that t is the parameter in the parametrization (j)(t ), and a is a 
new auxiliary variable; a similar construction was used in [5], see Definition 
3.2 therein, to analyze the rationality of the offset. 

We will restrict to the case when At is irreducible. Notice that At is 
reducible iff U 2 (t) +V 2 (t) is a perfect square, i.e. iff U 2 (t) +V 2 (t) is a perfect 
squar^J in which case OdiC) has two rational components. However this case 
admits an easier solution, as we observed in Section [271] Furthermore, as we 


5 In order to check whether or not U 2 {t) + V 2 {t) is a perfect square, we observe that 
non-negative elements of the ground held, i.e. non-negative real numbers, can be regarded 
as perfect squares. 
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Figure 3: Main idea to find the self-intersections of the offset. 


also observed in Section 2.1, if M. is irreducible then Od(C ) is irreducible, 
and therefore O d (C) is also simple. 

Now O d {C) can be seen as the image of A4 under the following rational 
transformation of the plane, denoted by <p(t,a) (see Figure [3j) : 


, , X(t) V(t) 

X = x(t, a) = 777777 + d~ 


W(t ) 


a 


y = y(t,a) = 


y(t) ,m 

W(t) a 


( 10 ) 


Since under our hypotheses O d {C ) is simple, for almost all points (x 0 ,y 0 ) G 
O d (C) there exists a unique ( t 0 ,a 0 ) such that (x 0 ,y 0 ) = tp(t 0 ,a 0 ). Therefore, 
the inverse ip^^x^y) exists for almost all points of O d (C), and is rational 
as well. Hence if 0d(C) is simple then tp defines a birational transformation 
between Ai and O d (C). 

Let us write 

V ) = y), a -1 (£, y)) , (11) 

where t _1 (a;,|/) and a~ l (x,y) are rational functions. Furthermore, let 


t \x,y) 


M x >y) 

B(x, y)' 


One can observe that if S G O d (C) is a self-intersection not generated by 
t — oo, then there are at least two different values ti,t 2 of t, and therefore 
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two different points Qi, Q 2 € Ai, such that <p(Qi) = ^{Qf) (see Figure[3]). So 
</9 _1 (S') cannot be defined, and therefore the denominator B(x, y ) of t _1 (x, y ) 
must vanish at S. Hence, the self-intersections of Od(C) are among the points 
of Od(C ) where B(x,y ) = 0. Since OdiC) is “parametrized” by tp(t,a), one 
does not need to know or make use of the implicit equation of Od(C) in order 
to find the points of OdiC ) where B(x,y ) = 0. Instead, one imposes 


B(x(t,a),y(t,a )) = 0. 


( 12 ) 


By repeatedly using that a 2 is a polynomial in the variable t, the numerator 
of (12) leads to a polynomial equation in t. The set of real roots of the 
polynomial provides a list of real t-values generating the real self-intersections 
of Od(C); notice that since OdiC ) is reduced, this list is finite. We will see 
later, in Section |322 , that in fact this list contains all the real non-isolated 
singularities of Od(C), not only the real self-intersections. 

In order to compute 0 1 (x, y), let us observe the following. For a particular 
(x 0 ,y 0 ) G Od(C), t 0 = t~ l (xo,yo) should be the unique root of 


gcd{P(x 0 ,y 0 ,t), Q(x 0 ,y 0 ,t)), 

for the polynomials P, Q introduced in Equation (|9]). Thus, in order to deter¬ 
mine the function t _1 (a;,|/), one can compute the gcd of P(x,y,t), Q(x,y,t ) 
for a generic point ( x,y ), considering P(x,y,t), Q(x,y,t ) as elements of 
M[ic,|/][t] (i.e. as polynomials in the variable t whose coefficients are real 
polynomials in x,y), with the additional condition F(x,y) = 0; recall that 
F(x, y) represents the implicit equation of the offset. More formally, one sees 
P, Q as elements of M(C ) c /(C))[t], where M((9 c /(C)) is the field of real rational 
functions of Od(C) (see [36j). Since OdiC ) is assumed to be irreducible then 
R(C?d(C))[i] is a Euclidean domain. Therefore 

Q(x,y,t) = gcd (P,Q) 

R(0 d (C))[t] 

is well defined and can be computed by means of the Euclidean algorithm. 
In [T7], Q{x,y,t) is introduced in a similar way, and is computed using the 
Euclidean algorithm. 


3.2. Computation of Q(x,y,t) and t 1 (x,y) via subresultants. 

In order to compute Q(x,y,t) by means of the Euclidean algorithm, we 
must perform several divisions between elements of M(C>d(C))[f]. Further¬ 
more, we need to check, at each step, if the remainder is zero in M(C^(C))[t]. 
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This implies either computing F(x, y ) first and then checking if F(x, y ) di¬ 
vides each remainder, which can be extremely costly, or substituting x = 
x(t, a), y = y(t , a) in each remainder and then checking if the result is divis¬ 
ible by a 2 - (U 2 (t) + V 2 (t)). 

In this section we will present an alternative, faster method based on 
subresultants. For this purpose we need several previous results. On the one 
hand, these results allow us to prove that Q(x,y,t), as a polynomial in the 
variable t with coefficients in x, y, has just one root; furthermore, we will see 
that this root is simple, so that the degree of Q(x,y,t) in the variable t is 1. 
On the other hand, these results will be used later to show that our method 
provides all the real singularities of Od(C), not only the self-intersections. For 
the next lemma we recall the notation X(t) = y(t) = introduced 
in Subsection [2] 

Lemma 8. Let (a?o, 2/o) £ Od(C), (xo,yo) f P± oo, which is not a self¬ 
intersection of Od[C), be an affine point generated by a real value t 0 with 
(F(t 0 ),V(t 0 )) f (0,0). Then ^(x 0 ,y 0 ,t 0 ) = 0 if and only if (x 0 ,y 0 ) is a 
local singularity of Od(C). 

Proof. Let us assume that (xo, yo) belongs to the exterior offset; we can argue 
in a similar way if (xo, yo) belongs to the interior offset. Since (xo, yo) G Od(C) 
and (xo,yo) is generated by t 0 , then 

P(xo,y 0 ,t 0 ) = Q(x 0 ,y 0 ,t 0 ) = 0 


and therefore 


Furthermore, let 


P(xo,yo,t 0 ) = Q(x 0: y 0 , t 0 ) = 0. 


P*(x, y, t) = (x- X(t)) ■ X\t) + (y- y(t )) • y\t). (13) 


A direct calculation shows that 

p ^ x,Vlt) = wW) ' p{ - XlVlt) - 


(14) 


Furthermore, since (xo, yo) is affine W(to) f 0. Since additionally P(x o, yo, t 0 ) 
0, we get P*(xo,yo,to) = 0. Furthermore, by differentiating in (14) we have 


dP* —3 W'(t) 1 dP 

~df~ w 4 (t) ( x,y,t ^ + wyf)'~dt' 
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Hence, since P(x 0 , y 0 , t 0 ) = 0 arid W(t 0 ) ^ 0, ^-(x 0 ,y 0 ,t 0 ) = 0 iff |y(x 0 , y 0 , t 0 ) = 
0. So let us see that, under the considered hypotheses, ^-(x 0 , ?/o, A>) = 0 iff 
(xo,Vo) is a local singularity of O c i(C). 

By differentiating (13) with respect to t, and evaluating at (xo,yo,to), we 

(xo,yo,t 0 ) = (xo~X (t 0 ))-X" (t 0 )+(yo-y(to))-y" (to) — ((X' (to)) 2 + (3^ (to)) 2 ) ■ 


get 
dP * 


dt 

Thus ^§j~(xo,yo,to) = 0 is equivalent to 

(*o - X(t 0 )) ■ X"(to) + (yo - y(to)) • y"(to) = (X’(to)) 2 + (y'(t 0 )) 2 . (15) 


Since (U(t 0 ),V(to)) ^ (0,0) by hypothesis, (X'(to), y'(to)) ^ 0- Therefore 
the normal vector to C at the point (X(t 0 ), y(to)) is ±(3^(to), — X'(to)), where 
we consider the + sign if (x 0 , yo) belongs to the exterior offset, as it is our case, 
and the — sign if (x 0 ,y 0 ) belongs to the interior offset. Additionally, since 
Q(xo,yoXo) — 0, from the definition of Q(x,y,t) we get that the modulus of 
(xq — X(t 0 ),y 0 — y(t 0 )) is equal to d. Hence under our hypotheses, 


(x 0 - x(t 0 ),y 0 - y(t 0 )) = 


d 


V*' 2 (t 0 ) + y 2 (t 0 ) 


(y'(to), —x'(t 0 )). (i6) 


The right hand-side of (15) is equal to the dot product of (xq — X(t 0 ),yo 


y(to)) and (X"(t 0 ),y"(to)). Therefore, taking (16) into account, we get that 

dP* 
at 


^-(xoiVoyto) = 0 is equivalent to 


x’(to)y"(to) - x"(t 0 )y'(t 0 ) 


[x >2 (to) + y / 2 (t 0 )] 


3/2 


1 

d' 


(17) 


This equality can be written as k(t 0 ) = where k(to) is the curvature 
of C at the point (X(t 0 ), y(to))- But under our hypotheses, k(t 0 ) = —^ is 
equivalent to (xo, yo) being a local singularity of Od(C) (see page 163 of [16]). 

□ 


Let (x 0 , yo) e O d (C), and let G Xom (t) = gcd(P(x 0 ,yo,t),Q(x 0 ,yo,t)). 
The roots of G Xo ^ yo (t) are exactly the affine t -values generating (xo,yo) via 
<j) d (t). The behavior of G XOtVO (t) is analyzed in the following result, which 
follows from Lemma [8] 
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Proposition 9. Let {x 0 ,y 0 ) G Od{C), (x 0 ,y 0 ) ^ P ±00 , be generated by some 
real t-value t 0 satisfying (U(to),V(to)) ^ (0,0). If (xo,yo) is a singularity of 
O d (C) then deg (G XOtyo (t)) > 1. 


Proof. If (xo, 7/o) is a self-intersection of Od(C), since by hypothesis (:ro, yo) 7^ 
P±oo, it must be generated by at least two t-values; therefore, deg(G iTO!yo (t)) > 
1. Suppose now that {x 0 ,y 0 ) is a local singularity of Od(C), not a self¬ 
intersection. Since {x 0 ,y 0 ) is generated by t 0 , it follows that t 0 is a common 
root of P(xo,yo,t) and Q(xo,yo,t). Moreover, by Lemma [8] 


dP 

P{x 0 ,y 0 ,t 0 ) = -^{xo,yo,to) = 0 . 


(18) 


Following the same argument as in the proof of Lemma [8j we get that (18) 
is equivalent to 

dP* 

P*i x o, 7/o, to) = -Qj-{x 0 ,yo,t 0 ) = 0. 

Additionally, Q(x,y,t ) is the primitive part of the numerator of 


Q*(x, y, t) = (x - x(t)) 2 + (y - y(t)) 2 - d 2 , (19) 

and = 2P *; in fact, one can check that the numerator of Q*(x,y,t) is 
Q(x,y,t). Hence, 


dO * 

Q*(xo,yo,t 0 ) = -^j-(x 0 ,y 0 ,t 0 ) = 0 . ( 20 ) 

Furthermore, since W(t 0 ) ^ 0 the content of Q(x,y,t ) cannot vanish at 
t = t 0 , and therefore ( |20| implies Q{xo,yo,to) = ^(xo,yo,to) = 0. As a 
consequence, to is also a root of both ^(xo,yo,t) and ^(xo,yo,t). Thus the 
multiplicity of t 0 as a root of G xom (t) is greater than 1. □ 

Remark 2. If(x 0 ,y 0 ) ^ P± oc is a real, non-isolated point ofOdiC), whenever 
deg(G xom (t)) > 1 then (x 0 ,y 0 ) must be a singularity of the curve hi defined 
by H(x,y) = Res t (P(x,y,t),Q(x,y,t)). Indeed, if G Xam (t) has just one 
multiple root to then since G XOiyo {t) is a real polynomial, t 0 € R; but then 
(xo,yo) must be a local singularity because of Lemma ^ If G X0iV0 (t) has 
different roots then ( x 0 , t / o ) is a self-intersection of'LL. By Proposition \23\ in 
Appendix III, the component of'LL defining Od(C ) is square-free, i.e. F(x,y ) 
(see page 7) is the polynomial of minimum degree defining Od(C). Therefore 
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if G X0t y 0 (t), with (x 0 ,y 0 ) G OdiC), has different roots then {x 0 ,y 0 ) is either 
a self-intersection of OdiC), or the intersection of OdiC) with some spurious 
factor (see Section 2.1). So z/deg(G X0i2/0 (t)) > 1 the point (xo,yo) is either a 
singularity of OdiC), or an intersection point between OdiC ) and a spurious 
factor of H. Hence, the converse of Proposition [5] is not true, in general. 

From Remark [2j we get that deg(Gb 0iyo (t)) = 1 for almost all points of 
( x o,yo) £ OdiC). Hence, the following result follows. 

Theorem 10. The degree of Q(x,y,t) in the t variable, is equal to 1. 

From Theorem [6] and Theorem 10, one has that 

Gix,y,t ) = Subresi(R,n, Q,m)(t), 


where V, Q represent the polynomials P , Q seen as elements of M((Td(C)), and 
n,m are the degrees of V, Q in M(C^(C)). However, from Lemma [l] we have 
n = deg t {Pix,y,t)) and m = deg t iQfx,y,f))] in other words, the leading 
coefficients of the polynomials Pix,y,t) and Q(x,y,t) with respect to t are 
not multiples of F(x,y), the implicit equation of OdiC). Therefore, by The¬ 
orem [7j in order to compute G{x,y,t) we can compute Subres i{P,n,Q 1 m) 
in the domain M[x, ?/][£], and then consider the coefficients of the resulting 
polynomial modulo F(x,y). After writing 


Subres! (P, n, Q, m){t) = sres^x, y) t + sr(x, y), 


we arrive at the following result. 

Theorem 11. The inverse mapping 0 satisfies 

.-if \ sr (x,y) 

t (x,y) = - 7 -r 

sresRx, y) 


( 21 ) 


for (x,y) G O d {C). 

Proof. Let Q(x,y,t) — B(x,y)t+A(x,y). Then t~ x {x,y) = — . However 

since G(x,y,t) = Subres i(P,n,Q,m)it) mod F(x,y), we have 


B(x,y) = sresi {x,y) mod F(x,y) 


and 


A{x,y) = sr fx,y) mod F(x,y). 


But then for (x,y) G OdiC), i.e. whenever Fix,y) = 0, we have (21). Oi 

Remark 3. The function y) is not necessarily defined at every point of 

OdiC). In fact, this is the crucial idea to find the self-intersections of OdiC). 
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4. Main result and complete algorithm. 

The following result, derived from the ideas and results in Section [3j 
shows how to compute the real singularities of Od(C) not coming from local 
singularities of C. 

Proposition 12. If (x 0 ,y 0 ) G O^iC.) is a non-isolated, real affine singularity 
of O d [C), generated by t 0 G M, with (U(t 0 ),V(t 0 )) f (0,0), then there exists 
a 0 G R, = U 2 (t 0 ) + V 2 (t 0 ), such that (x 0 ,y 0 ) = (x(t 0 ,a 0 ),y(t 0 ,a 0 )) and 

sresfxfo, a 0 ),y(t 0 , a 0 )) = 0. (22) 


Proof. Let (x 0 ,y 0 ) be a singularity of Od(C). Since ( x 0 ,y 0 ) is generated 
by a real value to, there exists «o such that (fo,«o) in M. and (xo,yo) = 
(x(t 0 ,ao),y(to,a 0 )). If (x 0 ,y 0 ) f P± oo, then we have deg(G X0! y 0 (t)) > 1 
by Proposition [9j On the other hand, by Lemma [l] we deduce that either 
deg(P(x 0 , yo, t)) = deg t (P(x,y,t)) or deg(Q(x 0l y 0 ,t)) = deg t (Q(x,y,t)). 
Furthermore, from Theorem [7] and Lemma [lj the subresultant chain of P. 0 
specializes well for x = Xo,y = yo (up to, perhaps, a constant). From The¬ 
orem |6j we can compute G XOtyo (t) as the first nonzero subresultant in this 
sequence. Since deg(Gj; 0i;yo (t)) > 1, we deduce that sresi(xo, yo) must vanish. 
If {x 0 ,yo) = P± oo, which means that P± oc is also generated by at least one 
real t, then sres^cco, Vo) = 0 by Lemma [I] and Definition [2j □ 


Therefore, if C does not have any local singularities, Proposition [12] pro¬ 
vides a method to find all the singularities of the offset. However if C has 
local singularities, then it remains to check that the singularities of OfC ) 
coming from local singularities of C also satisfy Equation (22). A proof of 
this fact is given in Appendix II. As a consequence, the following theorem, 
which is the main result of the paper, holds. 


Theorem 13. If (xo,yo) £ Od{C) is a non-isolated, real affne singularity of 
Od(C), generated by t 0 G M, then there exists «o G M, «o = ^ 2 (^o) + V 2 (fo), 
such that (x 0 ,y 0 ) = (x(t 0 ,oio),y(to,a 0 )) and sresi(a;(to, «o), 2/(^o, «o)) = 0. 

From Theorem [13] the real t values giving rise to the real, non-isolated sin¬ 
gularities of Od(C) are among the solutions of the system: 

sresi (x(t, a),y(t, a)) = 0, a 2 = U 2 (t ) + V 2 {t). 
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Substituting a 2 = U 2 (t ) + V 2 (t ) in sres^a^f, a), r/(t, a)) = 0 leads to a 
polynomial equation of the type 


fi(t)a + 77i(t) = 0. (23) 

Now squaring and using again that ct 2 = f/ 2 (t) + H 2 (f), we arrive at a uni¬ 
variate polynomial equation 

u(t):=m(&{t) + V*(t))-Tj*(t) = 0. (24) 


Note that the polynomial u)(t) cannot be identically zero because of Theorem 
10 Let u*(t) be the square-free part of u)(t), and let ui(t) = gcd ^( t f w(t)) • 
Since u(t) is not identically zero, the set B of real roots of cu(t) is finite. 
Hence, we deduce the following algorithm OffsetSing to compute a finite 
set containing the Lvalues generating the real, non-isolated singularities (self¬ 
intersections and local singularities) of O d {C). The algorithm requests O d (C) 
to be irreducible, which happens iff U 2 {t ) + V 2 {t) is a perfect square. 


Algorithm OffsetSing 

Require: A proper parametrization 0(f) of a planar curve C, and an offset¬ 
ting distance d > 0 such that O d (C) is irreducible. 

Ensure: A finite set B containing the Lvalues generating the real, non¬ 
isolated singularities of O d (C). 

1 : Let P(x,y,t), Q(x,y,t ) be the primitive parts of Q, ([5]), respectively. 

2 : Let Subresi (P,Q)(x,y,t) = sresi(x, y) t + sr(x, y). 

3: Substitute x = x(t, a), y = y(t, a) in sresi(a:, y). 

4: Find the real solutions for f of sresi(</?(t, a)) = 0, where a 2 = U 2 (t ) + 
V 2 {t). For this purpose: 

a. Compute the polynomial £i(t) a+rji(t) by substituting a 2 = f/ 2 (f) + 
V 2 {t ) in sresi(x(f, a), y(t, a)), and keeping the numerator. 

b. Compute the polynomial £ 2 (t) ( U 2 (t ) + V 2 (t )) — ?/ 2 (t), and u(t). 

c. Find the real roots of co(t). 

5: Return the list B of the real roots of u(t). 


Notice that squaring in £i(t) ■ a = which results from Equation 

(23), does not introduce fake solutions, because we are interested both in 
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the case a = J U 2 {t ) + V 2 {t ) and a = —-\Ju 2 (t) + V 2 {f)\ in fact, both 
correspond to offset points, one belonging to the exterior offset and the other 
one belonging to the interior offset. The only superfluous t values that can 
appear when computing the real roots of uj(t) are a consequence of Remark [2] 
These superfluous values appear when H[x,y ) = Res t (P(x,y,t),Q(x,y,t)) 
has extraneous factors that intersect Od(C), and correspond to the Rvalues 
generating these intersection points. Note also that if <f{t ) is polynomial, 
then there are no extraneous factors, and therefore no superfluous Rvalues 
appear. 

In order to identify superfluous t values, we first recall that the conditions 
for the appearance of extraneous factors, in the case of a non-polynomial <p(t), 


are described in Subsection 2T If these conditions hold, then the extraneous 
factors can be computed |15j , and any superfluous t value must give rise to a 
point on one of these extraneous factors. So in order to check if a real root t t 
of u>(t) is superfluous, one must first check if t t gives rise to a point on some 
extraneous factor. In the affirmative case, ti could still be non-superfluous 
whenever it generates a singularity (local or self-intersection) of O d (C). In 
order to detect this last situation, one can directly test: (a) if there exists 
another t 3 ^ ti, tj a real root of co(t), such that <f> d {ti) — 4>difj) (in which case 
tj and ti give rise to a self-intersection of G d (C)); (b) if lin \ t ^ t .(j)' d {f) = 0, in 
which case t — ti is a local singularity of O d (C). In practice, and unless ti is 
rational, (a) or (b) can be tested only up to a certain tolerance. 

Example 1. Let C be the cardioid, parametrized by 


X(t) = 


-1024 1 3 


256R + 32R + 1’ 


y(t) = 


-20481 4 + 128 1 2 
256R + 32R + 1' 


Let us now determine the singularities of the offset O d (C), for d — 1. We 
first observe that gcd(X(t), Y(t), W(t)) = l, and that the parametrization is 
proper. Moreover the polynomials U and V are 


U(t ) = 1024 1 2 (16 1 2 - 3) (161 2 + 1) , 

V(t) = -256 1 (48 1 2 - 1) (161 2 + 1) . 

The sum U 2 [t) + V 2 {t) = 65536 1 2 (16 1 2 + l) 5 is not a perfect square, and 
therefore O d (C) is irreducible for d — 1. Additionally, gcd (U(t),V(t)) = 
256 1 (16 1 2 + 1); hence, we define U(t ) = 4t (16 1 2 — 3) and V{t) = —48t 2 + l. 
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Now, the polynomials P(x,y,t) and Q(x,y,t) are 


P(x, y, t) = 64 xt 3 — (128 + 48 y) t 2 — 12 tx + y, 

Q(x, y, t ) = 256 ( x 2 + y 2 + 16 y + 63) t 4 +2048 xt 3 +32 (x 2 + y 2 - 8 y - l) f 2 +(x 2 +y 2 -l). 
J7ie principal subresultant sresi (x,y) is equal to 

sresi (x,y) = 1764 x — 5x 7 + 218 x 5 + 903 x 3 — 76x 5 y — 15x 5 y 2 

— 152 x 3 y 3 — 15 x 3 y 4 — 76 y 5 x — 5 y 6 x + 2596 yx 3 + 2020 xy 3 
+68 x 3 y 2 — 150 y 4 x + 10024 yx + 8823 xy 2 


After running the algorithm OffsetSing, the degree in t of the polynomial 
uj(t) is equal to 29. By computing its square free part and dividing out the 
common factors with W (t ) = 2561 4 + 32 f 2 + 1, we get 


(jj{t) = (t 4 + 


113 

9800 


r + 


12544 


3952 


3952 


t 4 + 


63232 


1011712 


Finally, by approximating the real roots of uj(t), we get B = {£i, t- 2 , £3, £4} 
where 


h = —0.04772, t 2 = 0.04772, t 3 = -0.08699, t 4 = 0.08699. 


The offset has a self-intersection at the points generated by i 3 and and two 
local singularities at the points generated by t\ and t 2 (see Figure^. Notice 
that C has a cusp, i.e. a local singularity, at the point (0,0), generated 
by t = 0. Ftowever, t = 0 does not belong to B. This is certainly not 
contradictory with Theorem 13. because one can check that the point (0,0) 
does not generate any singularity of Od(C). 


4-1. Growing of coefficients and degrees, and complexity of the algorithm. 

In this subsection we will analyze the growing of coefficients and de¬ 
grees in Algorithm OffsetSing, as well as the bit complexity of Algorithm 
OffsetSing. Notice that the bit complexity takes into account the growing 
of coefficients in the algorithm. We will use the standard Big Oh notation 
O , usually employed in complexity analysis, as well as the Soft Oh notation 
O , where logarithmic factors are ignored. Furthermore, in order to simplify 
the analysis we will assume that the parametrization has integer coefficients. 
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Figure 4: Cardioid curve 


We need to fix some notation: we denote by k the maximum of the degrees 
of X(t), Y(t ), W(t), and we denote by r the maximum of the bitsizes of the 
coefficients of X(t),Y(t),W(t). Recall that if the bitsize of the coefficients 
of a polynomial is r, then the coefficients of the polynomial are bounded by 
2 T . Furthermore, in our analysis we will assume that r >> log 2 (/c), so that 
log 2 (h) can be neglected when compared to r. This assumption corresponds 
to the case of highest computational cost. Notice that: 

(a) If we multiply two polynomials of degree k with coefficients bounded 
by 2 Tl and 2 T2 , the coefficients of the resulting polynomial are bounded 
by k2 T1+T2 . Hence, the bitsize of the product is bounded by T\ + r 2 + 
log 2 (fc). Assuming that Ti,72 >> log 2 (/c) the bitsize of the product 
of two polynomials, using the notation fixed at the beginning of the 
subsection, is 0{ji + r 2 ). If Ti = r 2 = r then we get O(r). 

(b) If we add r polynomials of degree k with coefficients bounded by 2 r , 
the coefficients of the sum are bounded by r2 r , so the bitsize of the 
coefficients of the sum is bounded by r + log 2 r. 

(c) By repeatedly applying (a), multiplying s polynomials of degree k with 
coefficients of bitsize r yields a polynomial where the bitsize of the 
coefficients is O(st). 

We also recall the following results on the complexity of some basic algo¬ 
rithms. We acknowledge here the help of Michael Sagraloff for pointing out 
some references and hints about the bit complexity of basic operations. 
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• The bit complexity of multiplying two univariate polynomials of degree 
n with coefficients of bitsize // is (D(ng) (see Corollary 8.27 of [3i]). 

• By repeatedly squaring and taking into account the complexity of mul¬ 
tiplying univariate polynomials, the computation of the p-th power of 
a polynomial f{t ) of degree n with coefficients of bitsize p can be done 
in 0{npp 2 ) time. 

• In order to multiply two polynomials F\ (t, s ) = sfi ( t ) + g\ ( t ) and 
F 2 (t, s ) = sf 2 {t) + g 2 (t), we can evaluate iq, F 2 at three random values 
of s, multiply the corresponding univariate polynomials, and then re¬ 
cover the value of F x (t, s ) ■ F 2 (t, s ) by interpolation. Hence if fi(t), p*(t) 
have degrees bounded by n and coefficients of bitsize p, then the cost is 
dominated by the cost of univariate polynomial multiplication, 0{np). 
In order to compute ( sfi(t ) + gi(t)) p we can apply a similar strategy. 
Since we need p + 1 evaluations of s, the bit complexity of computing 
Fi(t,s) ■ F 2 (t,s ) is the result of multiplying p + 1 times the bit com¬ 
plexity of computing the power of a univariate polynomial; hence we 
get 0(npp 3 ). 

• In order to compute the product of two bivariate polynomials F(t, s ) 
and G(t, s ) of degree n and bitsize p, we need to evaluate s at 2n + 1 
points si,..., S 2 n.+i, compute the products F(t, Sj) • G(t, s*), and then 
recover F(t, s ) ■ G(t, s ) by interpolation. Now F(t, Si ) and G(t, s*) are 
univariate polynomials of degree n and bitsize 0(n+p). Hence comput¬ 
ing F(t, ■ G(t, Si) has a cost of 0(n 2 + np). We need to carry out this 
process 2n + l times, so the total cost is 0(n 3 + n 2 /j,). The interpolation 
part is dominated by the total cost of the multiplications. Therefore 
the bit complexity of computing F(t,s) ■ G(t,s) is 0(n 3 + n 2 p); the 
coefficients of F(t, s) ■ G(t, s) have bitsize O(p). 

Now let us analyze each step of Algorithm OffsetSing. 

Step 1. The degrees of U(t), V(t) are bounded by 2k — 1, and the bitsizes of 
U(t), V{t) are, according to the above observations, O(r). Hence the degrees 
of P(x, y , t) and Q(x, p, t) are bounded by 3 k — 1, and their bitsizes are also 
0{t). In order to compute P(x,y,t) and Q(x,y,t) we need to remove the 
t-contents of P(x,y,t) and Q(x,y,t). By Lemma 11 in [23], a divisor of a 
polynomial in Zfaq ,... ,X(\ of degree N with coefficients of bitsize bounded 
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by /i has coefficients of bitsize 0(N + g). Since P(x,y,t ) and Q(x,y,t ) 
are divisors of P(x,y,t ) and Q(x,y,t), in our case N = 3k — 1 and fi = 
0 (t ); hence, P(x,y,t ) and Q(x,y,t ) have coefficients of bitsize C)(fc + r). 
Furthermore, the degrees of P(x,y,t ) and Q(x,y,t ) are bounded by — 1. 
So we finish this step with polynomials P, Q of degrees d?(fc), and bitsizes 
d?(Jfc + T). 

The bit complexity of step 1 is dominated by the computation of the 
/-contents of P, Q. From (jb]), Q, (J8|, the /-contents of P, Q are products of 
gcds of univariate polynomials of degree 0(k) and bitsize 0{r). Hence the 
bit complexity of step 1 is 0(k 2 T ) (see Section 11.2 in [31]). 

Step 2. Let us analyze the computation of Subresi (P,Q)(x,y,t). According 
to PS], if f,g G rf], deg t (/) = p, deg t (^) = q, deg y .(f) < e u 

deg (g) < 6i, e = ei • • • e^, and the bitsizes of /, g are bounded by /i, then: 

(i) The degree in y t of each subresultant of /, g with respect to t is bounded 
by (p + q)e.i. 

(ii) The complexity of the computation of any element of the subresultant 
sequence is 0(q(p + g) £+1 e/i). 

(iii) The total degree in yi ,..., ye of each subresultant is bounded by (p + 

q) Eli 

(iv) The bitsize of the coefficients of the subresultants is bounded by (p + 
q)p. 

In our case, / := P, g := Q, i = 2, y x := x, y 2 ■= y, p = q = 0(k). 
Furthermore, since deg a .(P) = deg ?/ (P) = 1, deg I ,(Q) = deg y (Q) = 2, then 
e\ = 2, e 2 = 2, and therefore e = 4. Also, p = 0{k + r). Hence, 

(i) The degrees in x or y of each subresultant of P, Q with respect to t is 
O(k). 

(ii) The complexity of the computation of any element of the subresultant 
sequence is 0(k 5 + k A r). 

(iii) The total degree in x,y of each subresultant is 0(k). 

(iv) The bitsize of the coefficients of the subresultants is 0(k 2 + kr). 
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Therefore the bit complexity of step 2 is (D(k 5 + k 4 r). 

Step 3. Let Subresi(P, Q)(t) = sresi(x,y)t + sr(x,y). According to the 
previous analysis, sresi (a;, y) is a bivariate polynomial of degree N where 
N = 0(k). Therefore the number of terms of sresi (x,i/) is bounded by 
( JV + 2 ) = ^(^ 2 )- Now we need to substitute x := x(t,a ) and y := y(t,a) in 
sresi(a;, y). Notice that 


Writing 


x(t, a) = 
y(t,a) = 


m 

w (t) 
Y(t ) 
W(t) 


+ d 

a 

a 


V(t) 


aX (t) + dV(t)W{t) 
aW(t) 

aY(t) - dU(t)W(t) 
aW(t) 


sresi (x,y)= ^ a i:j x l y\ 

0<i,j<N 
0 < i + j < N 


we have that sresi (x(t,cn),y(t,a)) is equal to 


(aX(t) + dV(t)W(t)Y ■ ( aY{t ) - dU(t)W(t)) j ■ [aW(t)} N ~^ +j) 

^ aij a N W N {t) ' 

0<i,j<N 

0 < i + j < N 

(25) 

The total degrees, as polynomials in a,t of aX(t) + dV(t)W(t), aY(t) — 
dU(t)W(t ) are bounded by 3 k — 1, and the total degree of aW(t) is bounded 
by k + 1. Hence the total degree of each term of the numerator of (25) is 
bounded by 


(3 k - l)(i + j) + (k + 1)(N - i - j ) = (2k - 2 )(i + j) + (k + 1)N, (26) 


where 0 < i + j < N. Thus, the total degree of (25) is bounded by the result 
of replacing i + j — N in ( J26| ). Since N = 0(k ) we get that the total degree 
of the numerator of (25) is 0(k 2 ). Notice also that the degree in a of the 
numerator of (j25j) is bounded by N = 0(k). 

Since U(t),V(t) are factors of U(t),V(t) and the bitsize of the coefficients 
of U(t),V ( t ) is 0(t) and their degrees are 0(k ), then the bitsize of U(t),V ( t ) 
is 0(k + r). Then the bitsize of the coefficients of (aX(t) + dV(t)W(t ))* 


26 












is 0(i(k + r)), the bitsize of the coefficients of (aY(t) — dU(t)W(t)y is 
0(j(k + r)), and the bitsize of the coefficients of [aW is bounded 

by 0{(N — i — j)r). For each i,j we have i + j + [N — (i + j)] = N, and 
therefore the bitsize of the coefficients of (aX(t) + dV{t)W{t)) 1 ■ ( aY(t ) — 


dU(t)W{t)y ■ [aW{t)} N ~^ + ^ is 0(N(k + r)) = 0{k 2 + kr). Since the bitsize 
of a t j is 0{k 2 + kr), the multiplication by yields coefficients of bitsize 
0({k 2 + kr) 2 ). We may need to add 0{k 2 ) coefficients of this bitsize, but 
neglecting logarithmic terms we again get bitsize 0((k 2 + kr) 2 ). 

Let us compute now the complexity of step 3. We observe that each 
(aX(t) + dV(t)W(t ))* and ( aY(t) — dU(t)W(t)) j in (26) has degree 0{k 2 ) 
and bitsize 0{k 2 + kr). Therefore from the results at the beginning of the 
subsection, the complexity of computing each multiplication 


(. aX(t ) + dV(t)W(t)Y ■ ( aY(t ) - dU(t)W(t )) j • [aW(t)} N ~ {i+j) 


in (26) is dominated by O ((k 2 ) 3 + (k 2 ) 2 (k 2 + kr)) = 0(k 6 + k 5 r). Since 
sresi(a;, y) has 0(k 2 ) terms, we need to perform 0(k 2 ) multiplications of 
this kind. Hence we get a complexity 0{k 2 ■ ( k 6 + k 5 r)) = 0(k 8 + k 7 r) for 
step 3. 


Step 4■ Let a(t,a ) be the numerator of sresi (x(t,a),y(t,a),t), i.e. 

a(t, a) = cza r(t)a N + aAr_i(t)a Ar_1 + • • • + a 0 (t). 


Recall that the total degree of a(t, a) is 0(k 2 ), and the degree in a of a(t, a) is 
0(k). Now we need to substitute a 2 := b(t) = U 2 (t) + V 2 (t) into a(t, a). The 
degree of b(t) is 0(k), and the bitsize of the coefficients is 0(k+r). Hence the 
computation of each power ( U 2 (t ) + V 2 {t)y has complexity 0((k 2 + kr)] 2 ). 
Since j < k, this complexity is dominated by (9(fc 4 + k 3 r). Furthermore, the 
bitsize of the coefficients of each (U 2 (t) + H 2 (t)) J is dominated by 0[k 2 + kr), 
and the degree in t of ( U 2 {t) + V 2 {t)y is dominated by 0(k 2 ). The degree in t 
of each a 3 (t) is Oik 2 ), and the bitsize of its coefficients is 0((fc 2 +fcr) 2 ). Hence 
the complexity of each multiplication aj(t)(f/ 2 (t) + V 2 (t)y* is dominated by 
0{k 2 ■ (k 2 + kr) 2 ) = 0{k A {k + t) 2 ), and gives rise to coefficients of bitsize 
0((k 2 + kr) 2 ). We need to compute 0{k) products of this kind, so we get a 
total complexity of 0(k ■ k A {k + r) 2 ) = 0(k 5 (k + r) 2 ). The degree in t that 
we get is 0{k 2 ). 

Finally we arrive at equation (23) and then, after squaring, we compute 
the polynomial &{t) in (24), and then oj{t). These operations do not modify 
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the total complexity of step 4. Furthermore, the degree of oj{t) is 0[k 2 ), and 
the bitsize of the coefficients of oj(t) is 0((k 2 + kr) 2 ). According to Theorem 
5 in pa, the cost of isolating the roots of uo(t) is 0(k 6 (k + r) 2 ). Therefore, 
the overall bit complexity of step 4 is 0(k 6 (k + t) 2 ). 

The computational cost of checking whether or not A4 is irreducible 
is very low, and does not modify the overall complexity of the algorithm. 
Adding the complexities for step 1, step 2, step 3 and step 4, we achieve 
an overall bit complexity of 0(k 6 (k + r) 2 ) for Algorithm OffsetSing. We 
summarize the previous reasonings in the following result. 

Theorem 14. The overall bit complexity of Algorithm OffsetSing is bounded 
by 0(k 6 (k + t) 2 ). The growing of the coefficients in the algorithm is domi¬ 
nated by 0((k 2 -\-kr) 2 ), and the growing of the degrees is dominated by 0{k 2 ). 
Furthermore, the polynomial u{t) computed by Algorithm OffsetSing has 
degree 0(k 2 ) and coefficients of bitsize dominated by 0({k 2 + kr) 2 ). 

4-2. Experimental Results 

We have implemented our algorithm in the computer algebra system 
Maple 16. Since the algorithm requires to compute subresultants, we tested 
several algorithms for carrying out this operation (see P); in practice, the 
best method seems to be the direct computation of the subresultant chain 
using the command RegularChains [ChainTools] [SubresultantChain] of 
Maple 16. Some of the rational curves we have tested, denoted by Cj, 
i G {1,..., 8}, are enlisted below: 


_ t +1 3 _t~t 3 

~ 1 + t 4 ' y - 1 + t 4 

_ -7f 4 + 288f 2 + 256 _ -80f 3 + 256f 

“ t 4 + 32 1 2 + 256 l V ~ t 4 + 32f 2 + 256 


C 3 


C 4 


C 5 


C 6 


c 7 


C 8 


18 f 4 + 21 f 3 — 7f — 2 36f 4 + 84t 3 + 73f 2 + 28f + 4 

“ 18 f 4 + 48 f 3 + 64 f 2 + 40 f + 9 ’ ^ _ 181 4 + 48t 3 + 64t 2 + 401 + 9 

x = l- (l + t 2 )' 1 ,y = t- 

-t 4 - 6 1 2 + 3 0 f 3 

a; =- k—, y = 8-o 

(1 + t 2 ) 2 (1 + t 2 ) 2 

1 - 3f 2 _ (1 - 3t 2 )t 

X “ (1 + P) 2,V ~ (1 + f 2 ) 2 

t 2 — 3 t(—t 2 + 3) 

1 + t 2 ’ y 1 + t 2 

_ 87 - 7 1 4 + 22t 3 - 55t 2 - 94 1 _ -82 - At 4 - 83 1 3 - 10f 2 + 62 1 

-73 - 56t 4 - 62f 2 + 97t ’ V ~ -73 - 56t 4 - 62f 2 + Wt 


The results of our experiments on the curves C* are summarized in Ta¬ 
ble 1: here, the variable Time shows the computing time (in seconds); d is 
the offsetting distance; n p is the number of real /-values computed by the 
algorithm; S t is the degree of ou(t), r is the bitsize of u(t), and 5 t (P),S t (Q ) 
are the degrees in t of P(x,y,t ) and Q(x,y,t). Additionally, we include the 
degree deg (F(x,y)) of the implicit equation of the offset, in order to give an 
idea of the size of the problem in each case. All the computations have been 
carried out with the Computer Algebra System Maple 16 in an iMac with 
an Intel Core i3 processor with speed revving up to 3.06GHz. 
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Example 

Time 

d 

n p 


T 

8t(P)MQ) 

S(F(x,y)) 

Ci 

0.140 

1 

10 

30 

22 

6, 4 

12 

C 2 

0.122 

1 

9 

21 

41 

4,4 

10 

C 3 

6.167 

1 

26 

222 

510 

10, 8 

20 

C 4 

0.051 

1 

4 

22 

16 

4,4 

8 

C 5 

0.068 

1 

8 

20 

23 

3, 6 

8 

C 5 

0.092 

0.3 

12 

22 

53 

3, 6 

8 

C 6 

0.244 

1 

21 

81 

84 

6, 6 

14 

C 7 

0.074 

1 

9 

29 

26 

5,4 

10 

C 8 

6.948 

0.5 

12 

228 

927 

10, 8 

20 


Table 1: Examples. 

The pictures corresponding to the examples in Table 1 can be found in 
Figure[5] from left to right, the curves Ci, C 2 , C 3 appear in the first row, C 4 , 
C 5 with d — 1, C 5 with d = 0.3 are shown in the second row, and C 6 , C 7 , C 8 
are shown in the third row. In all the cases, with the exception of C 7 , we get 
exactly the set of t -values generating the singularities of the offset, without 
any extra values. In the case of C 7 , we get a superfluous value, namely 
t = 0. This t-value generates the offset points (—4,0) and (—2,0). None 
of these points are singularities of the offset; however, the point (—2, 0) is a 
singularity of "H = Rest(P, Q ), which in this case has extraneous components, 
and corresponds to the intersection of one of these spurious components with 
the offset to C 7 . 

Notice that in general n p is bigger than the number of singularities, since 
several Rvalues give rise to the same self-intersection. For instance, in the 
case of C ,5 with d = 1, the algorithm provides 12 t- values. However, the curve 
has 9 singularities; 6 of these singularities are local singularities, each one 
generated by a different Rvalue, and three singularities are self-intersections, 
each one generated by two different Rvalues. Therefore, we get = 6 • 1 + 
3 • 2 = 12. 

4-3. Comparison with other approaches 

We have also implemented in the same computer algebra system the algo¬ 
rithm described in [15] , in order to compare timings. We provide the details 
of the experiments carried out with several polynomial curves, denoted as 
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Cj, i G {9,..., 13}, enlisted below: 


C 9 : x = 87 f 5 + 44 f 4 + 29t 3 + 98 f 2 — 23 f + 10, 
y = —61f 5 — 8 f 4 — 29 f 3 + 95 f 2 + 11 f — 49. 

C 10 : x = l + f 4 + 2f 3 + 9t 2 -3f, 
y = —4 + f 3 + 5t 2 + t 

C n : x = 95 f 5 + 11 f 4 — 49f 3 — 47f 2 + 401 — 81, 
y = 981 5 - 23t 4 + 101 3 - 61 1 2 - 8 1 - 29. 

C12 : x — t, 

y = t\ 

Cis : x= ^ (l-t) 6 + ^ (l-t) 5 t-15 (l-f) 4 f 2 + 15 (l-f)t 5 +-^-f 6 , 

10 5 5 10 

y = (l-tf+ 2 ^(l-t) 5 t + 9 (1 — t) 4 1 2 + 9 (l-t) 2 t 4 + ^ (l-t)t 5 + t 6 . 

0 0 


The following table, Table 2, includes the parameters in Table 1 pins Time 2 , 
which is the computing time of the algorithm in [15]. 


Ex. 

Time 

Time2 

d 

n p 

s t 

T 

St(P)MQ) 

S(F(x,y)) 

c 9 

0.830 

60.286 

2 

8 

200 

660 

9,10 

18 

c 10 

0.247 

4.527 

1 

4 

108 

177 

7,8 

14 

C n 

1.288 

60.443 

5/3 

4 

200 

696 

9, 10 

18 

c 12 

0.107 

0.667 

0.8 

8 

84 

99 

7,8 

14 

c 13 

4.672 

475.200 

0.05 

4 

320 

948 

11, 12 

22 

c 13 

6.779 

555.871 

~ 0.03141 

4 

320 

1547 

11, 12 

22 


Table 2: Comparative examples 
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The timings in Table 2 show that our algorithm is clearly much faster 
than the algorithm in [15]. The pictures corresponding to the examples in 
Table 2 can be found in Figure [6j from left to right, the curves Cg, Cio, Cn 
appear in the first row, Ci 2 , Ci 3 with d = 0.05 and Ci 3 with d =~ 0.3141 
are shown in the second row. The curves Ci 2 and C 13 are also considered in 
[21] : C 12 is the superbola y = t A and C i3 is the bottle-shaped Bezier curve 
whose control points are given by 

(-0.1,1), (0.3, 0.7), (-1, 0.6), (0,0), (1, 0.6), (-0.3, 0.7), (0.1,1). 

Note that the curve Ci 3 for d ~ 0.3141 is an example of a tacnode. Our 
algorithm computes four different real t-values, but with an appropriate tol¬ 
erance one observes that only two of these values are regarded as different. 
In [2TJ, several timings are listed for C i2 and Ci 3 , depending on the imple¬ 
mentation method (double precision floating point arithmetic, or rounded 
interval arithmetic), and the tolerance used in the computations. Although 
our timings are better than those in EH. here the comparison is less clear, 
since the method in [[21] is implemented in C++ and the timings correspond 
to a graphics workstation running at 36 MHz, very different from our own 
machine; implementing in our own system the technique used in [21] would be 
really difficult, since a number of nontrivial algebraic and numerical strategies 
are involved. 

5. Conclusions. 

In this paper we have presented a novel method to compute the real affine, 
non-isolated singularities of the offset of a planar curve described by means 
of a proper rational parametrization (X(t),y(t)). The method is based on 
ideas in na. and requires the computation of the inverse of a mapping that 
relates the offset with an auxiliary, simpler curve A4. This curve A4 lives 
in the plane (t,a) and has the equation a 2 = U 2 {t ) + V 2 (t), where t is the 
original parameter of the curve and a is an auxiliary variable. The method 
is easy to describe and to implement, and provides a finite list containing 
the affine t- values generating the real, non-isolated singularities of the offset, 
which is useful, in particular, for trimming applications. The method can be 
generalized to other geometric constructions involving square roots studied 
in the CAGD literature, like for instance generalized offsets (0, M or 

conchoids m 1321 , [ 33 ]). 
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Figure 6: Polynomial examples 
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6. Appendix I: proof of Lemma [l| 

The goal of this appendix is proving Lemma [TJ In order to do it, we need 
to analyze the behavior of P,Q at the points of Od(C) generated by t = oo, 
if any. We introduce the notation poo = Km WJt))- 

Lemma 15. The leading coefficient of Q with respect to t vanishes if and 
only if Poo is an affine point. Furthermore, the only points where the leading 
coefficient of Q with respect to t vanishes are the points of the circle of radius 
d centered at poo. 


Proof. Note first that the leading coefficient of Q with respect to t is equal 
to the leading coefficient of Q up to a multiplication by a constant. Dividing 
([5]) by W 2 (t), we get 


f(x,y,t ) := x 


m V 

W(t)J 


+ y 


m \ 2 

W(t)J 


d 2 


One can see that deg(IF) > max{deg(X), deg(Y)} iff Poo is affine. Further¬ 
more, if deg (IF) > max{deg(X), deg(F)} then deg(Q) = 2 ■ deg(IF), and 
the leading coefficient of Q with respect to t is the product of the leading 
coefficient of W 2 {t), multiplied by lim f(x,y,t). But this coincides with the 

£—» 00 

equation of the circle of radius d centered at P^. 

Finally, if deg(IF) < max{deg(X), deg(F)}, in which case poo is not affine, 
then the leading coefficient of Q with respect to £ is a nonzero constant. □ 


Lemma 16. The leading coefficient of P with respect to t vanishes iff Poo is 
an affine point. Furthermore, the leading coefficient of P with respect to t 
vanishes over the line normal to C at Poo- 


Proof. First, note that the leading coefficient of P with respect to t is equal 
to the leading coefficient of P , up to a multiplication by a constant. Let 
r := max{deg(X), deg(F)}, s := max{deg(f/), deg(F)} and w := deg(IF). 
Furthermore, let c(W, r) and c(W, w ) denote the coefficient of t r and t w re¬ 
spectively in IF(£); similarly for X and Y. Also, let c(U,s ) denote the coef¬ 
ficient of t s in U{t)\ similarly for c(V , s). 

Now Poo is affine iff deg(IF) > r. Furthermore, if deg(IF) > r then the 
leading coefficient of P with respect to t, is: 


lcoeff (P) = c{U , s) ■ [c(W, w)-x- c(X, w)] + c(F, s) ■ [c(IF, w) ■ y - c(Y, w)]. 
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This expression can be written as the product of the constant c(W, w ) • 
\JC 2 (U, s ) + C 2 (V, s ) times the dot product of the following two vectors: 


Vl 


c(U, s ) 


c(V,s) 


V / c 2 (f/,s) + c 2 (^s) , y/ c 2 (U, s) + c 2 (V, s) 


and 

v 2 = (x, y) 
Furthermore, one can see that 


( c(X,w) c(y, w) \ 
\c{w, w y c{w,w ))' 


„ _ lini , m _ v(t) 

Vl \y/U 2 (t)+V 2 (tY y/lP{t) + V 2 (t) i 

which is the limit of the unitary vector tangent to C at p Additionally, 
denoting p^ we also have that v 2 = (x,y) - (loo,^), which 

represents the vector connecting p^ and a generic point (x,y). Therefore, 
V\ ■ v 2 — 0 is the equation of the line normal to C at p^. 

If deg (IF) < r, in which case p^ is not affine, one can easily see that 

lcoeff(P) = — (r — w)C(W, w)(C(X, r) 2 + C(Y, r) 2 ), 


which is always nonzero. 

Then we can finally prove Lemma [TJ 


□ 


Proof, (of Lemma [I]) From Lemma 15 and Lemma 16, we have that the 
leading coefficients of P(x,y,t), Q(x,y,t ) with respect to t simultaneously 
vanish at the points (xo,yo) E Od(C ) belonging, at the same time, to the 
circle of radius d centered at p^, and to the line normal to C at p^. But 
these are exactly the points of O^fC ) generated by p^. □ 


7. Appendix II: singularities of Od(C -) coming from local singular¬ 
ities of C. 

The goal of this appendix is to prove that the singularities of Od(C) 
coming from local singularities of C also satisfy Equation (22). In this sense, 
if (x 0 ,y 0 ) E Od(C ) is generated by t 0 6 M, with (U(t 0 ),V(t 0 )) = (0,0), and 
is a self-intersection, then sresi(;r(fo, Q?o), y(to, cto)) = 0. Observe that if 
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{x 0 ,y 0 ) ~f~ P± oo then deg(G a;0i y 0 (t)) > 1 and so, sresi(x(t 0 , a 0 ),y(t 0 , a 0 )) = 0; 
otherwise, (xo,yo) = P± oo, sresi(a;(to, «o), y(to, «o)) = 0 by Lemma [I] and 
Definition [3 

Hence, in the rest of the appendix we will focus on the case when (x 0 , yo) 
is not a self-intersection. 


Lemma 17. Let (xo,yo) E Od(C) be generated by to E M, and assume that 
(U(to),V(to)) = (0,0). Then the multiplicity of t 0 as a root of G Xo>yo (t) is 
higher than 1 if and only if 

U{t 0 )V'{t 0 ) - V(to)U\to) = 0. 


Proof. We observe first that since W(t 0 ) ^ 0, to is a r °ot of Q(xo,yo,t) 
of multiplicity k iff t 0 is also a root of Q(xo,yo,t) of the same multiplicity. 
Similarly, t Q is a root of P(x o, yo , t) of multiplicity £ iff to is also a root of the 
polynomial 


P(x 0 ,yo,t) 


PjpyP) 

gcd (U{t),V{t)) 


of the same multiplicity. Furthermore, one can easily see that since W (to) 7^ 
0, if t 0 is a root of P(x 0 ,yo,t ) of multiplicity £ then t Q is also a root of 
Q( x o,yo,t), and therefore of Q(x 0 ,y 0 ,t), of multiplicity £ + 1. So it suffices 
to show that the condition in the statement of the lemma holds if and only 
if the multiplicity of t 0 as a root of P(xo,yo,t) is higher than 1. In order to 
do this, notice that the derivative of 


P(x o, Vo, t)/W (' t ) = U(t) ( xq - + V(t) (vo - ^7^ 


is: 


„,u, l A '(() \ ( Y(t) \ U(t)U(t) + V(t)V(t) „„ 

- mo —- (27) 


Then t = t 0 , where W(to) ^ 0, is a root of P(x 0 ,y 0 ,t) of multiplicity higher 


than 1 if and only if (27) vanishes at t — t 0 ■ Since U(t 0 ) = V(t 0 ) = 0 by 


hypothesis, the evaluation of (27) at t — t 0 is equal to the dot product of the 
vectors: 

d=(Ufto),Vfto)),b= ( Xo -§^,y 0 - F(t0) 


W(t 0 


W(to 
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Since (x 0 ,y 0 ) *= 0^(0) is generated by t — t 0 , we have that b is parallel to 

V(t 0 ) U(t 0 ) 


U^t 0 ) + v^t 0 ) Jir>(t 0 ) + v*(to) 


□ 


Thus, (27) vanishes at t = t 0 if and only if 

U(t 0 )V'(t 0 ) - V(t 0 )U'(t 0 ) = 0. 

The following corollary follows from Lemma 0 

Corollary 18. Under the hypotheses of Lemma 11, U {to)V' {tf)—V {tfjU' {tf) = 
0 implies sresi(xo, yo) = 0. 


Proof. By Lemma [17] we have that deg (G X0M (t)) > 1. Then one argues as 
in the proof of Proposition [12] □ 

We need some previous work in order to see the geometrical meaning of 
Lemma 17 Let T : M 2 — y M 2 be an orthogonal change of coordinates 

T(x, y) = (ax + by + c x , —bx + ay + c 2 ), 
with a 2 + b 2 = 1. Let also 

Mt) = T(x(t),y(t)) = (x(t),y(t)), 
with X(t) = aX(t ) + by(t) + c\, y(t ) = — bX(t) + ay(t) + c 2 , and let 
T(U(t),V(t )) = (U(t),V(t)) and T(U(t),V(t)) = (U T (t),V T (t)). 
Hence, it is easy to see that 

r (t) = |W , y'(t) = 


W 2 (t)’ v 7 W 2 (t) 

Moreover, since 

gcd (U, V) = gcd(aU(t) + bV(t), -bU(t ) + aH(t)) = gcd(£7(t), V(t)), 
then 

C T (t) = U(t)/gcd(U(t),V(t)), V T (t) = V(t)/gcd(U(t),V(t)). 
Now we have the following instrumental lemma. 
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Lemma 19. Let 


i(t) = x'(t)y"(t)-x"(t)y'(t). 


and let 

i T (t) = r(t)y"(t)-x"(t)y'(t). 

The following statements are true: 

(1) m=£ T (t). 

(2) U T (t)Vf(t) - U' r {t)V T (t) = U{t)V'(t) - U'(t)V{t) 

Proof. (1) The equality can be verified by a direct computation. (2) Let 
v{t) = gcd (U{t),V{t)). Observe first that 

m vW (t) „ v(t) i4t)V(t) 

1 1 W 2 (t) W 2 (t) ’ y W W 2 (t) W 2 (t) ’ 

and therefore 

_ v\t) ■ (U(t)V'(t) - U'(t)V(t )) 
j W\t) 

Furthermore, 

f m V (t)v T {t) yw . ughw 

11 IV>(«) IV 2 (i) iv 2 (i) IV 2 («) ’ 


and hence 


It it) 


v\t)-(y T (t)vf(t)-u' T (t)v T (t)) 

W\t ) 


Finally, using the statement (1), the statement (2) follows. 


□ 


Before giving a geometric translation of Lemma [l7j we need an additional 
ingredient, namely the notion of place recalled in Subsection T3 So let S' G C 
be a local, real affine singularity of C, S = (p{to), to G M, and let us consider a 
coordinate system centered at S where the x-axis coincides with the tangent 
to V at S (see Figure [7]). Then S is the center of a place of C that can be 
written as 

V{h) = {h*rf q h* + •••), 


with p, q G N, p > 2, q > p. In [2], the problem of determining the places 
V± d (h) of Od(C ) generated by a given place V(h) of C under the offsetting 
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Figure 7: Coordinate system for Proposition 


20 


transformation is addressed. Furthermore, from Theorem 7 in [2], one can 
see that if S' is a local singularity of C, then it generates a local singularity 
of OdiC ) if and only if q — p > 2. Now let us see that Lemma 17 implies 
that S must generate a singularity of Od(C). This is done in the following 
proposition, where we use the previous notation and the ideas in Lemma [19] 


Proposition 20. Let ( x 0 , yo ) £ Od(C) be generated by t 0 e K, and assume 
that (U(t 0 ),V(t 0 )) = (0,0). Then (x 0 ,y 0 ) is a local singularity of OdiC) iff 
Uito)V'it 0 ) - V(top(t 0 ) = 0. 


Proof. (4=) Let S = 0(£o) G C, and let x' be the tangent line to C at S 
corresponding to t —> t 0 . Furthermore, let y' be the line perpendicular to x' at 
S', and let T be the orthogonal change of coordinates mapping the coordinate 
system {0]x,y} onto the coordinate system {S;x',y'} (see Figure [7]). Also, 
let (T(t), y(t )) = T(A(t), y(t )); notice that (T(t), y(t)) is a parametrization 
of C in the coordinate system {S'; x', y'}. Since 


(x\t),y t (t)) = ^-(u T (t)Mt)), 


we have that 

m(t) = 

X'(t) U T (t ) 
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Observe that the tangent to C at S is parallel to (Hr(to), Vr(^o)) f (0,0). 
Since this tangent is parallel to the .F-axis, Ur {to) f 0 and m(to ) is well- 
defined. Furthermore, 


m'(t) = 


Vf(t)U T (t) - ^(t)P T (t) 

om 


and by the condition (ii) and the statement (2) of Lemma 19, we have that 
m'(t 0 ) = 0. 

On the other hand, let us consider a place {x{h),y{h)) = {h p , j3 q h q + • • •) 
of C centered at S. Then we have that 


11 \ y'(h) q(3 q h q 1 + 

n(h) := = 


= h q - r + ■■ ■ 
P 


(28) 


Since ( x{h),y{h )) also parametrizes C around S, and is written in the same 
coordinate system as (X(t),y(t)), namely the coordinate system {S; x\ y'}, 
for every t sufficiently close to to we can find h such that 


y'{h) y'jt) 

x'{h) X'(tY 

Therefore n'(0) = m'(to ) = 0. However, n'(0) = 0 implies q — p > 2, which 
is the condition (see Theorem 7 in |2]) for S to generate a singular point in 
O d {C) . 

(=>) If (a: 0 , 2 / 0 ) is a local singularity then q — p> 2, and therefore n'{ 0) = 
m'(t 0 ) = 0; but then U(t 0 )V'(t 0 ) — V(t 0 )U'(t 0 ) =0. □ 

Finally, we can eventually prove the aimed result. 


Corollary 21. If {xo,yo) G Od(C) is a non-isolated, real affine singularity 
of Od{C), generated by t 0 G M, with {U{t 0 ),V{to)) = (0,0), then there exists 
«o G K, = U 2 {t 0 ) + V 2 {to), such that (xo,2/o) = (a; (to, «o), y(to, «o)) and 
sresi(x(t 0 , a 0 ),y(t 0 , «o)) = 0. 


Proof. If (a;o,2/o) is a self-intersection of Od{C), the result follows from the 
observations at the beginning of the appendix. If (x 0 , y 0 ) is a local singularity 
of Od(C), then the result follows from Corollary 18 and Proposition [20] □ 
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8. Appendix III: structure of Res t (P(x,y,t),Q(x,y,t)). 


From Section 2.1 we know that the polynomial 


H(x,y) = Res t (P(x,y,t),Q{x,y,t)) 


can be written as H(x,y ) = F(x,y) ■ G(x,y), where F(x,y ) = (f(x,y)) r , 
f(x, y) is an irreducible polynomial implicitly representing and G(x, y) 

is the product of all the extraneous factors. Our goal is to prove that, under 
the hypothesis that Od(C ) is simple and C is properly parametrized, we have 
r = 1, i.e. the component of H(x,y) corresponding to OfiC) has multiplicity 
1 . 

Let us consider the set A of the yos satisfying some of the following 
conditions. Here we use the notation P± oc = (x± 00 ,y± 00 ). 


(1) The intersection of the line y = yo with Od(C) contains some point also 
belonging to the curve G(x, y) = 0. 


(2) The leading coefficients of P(x,y,t ) and Q(x,y,t ) with respect to t 
identically vanish when y = y$. 


(3) yo = y± oo- 

(4) The line y = yo is tangent to Od{C). 

(5) The line y — y 0 contains either a local singularity of Od(C), or a point 
of Od(C ) simultaneously generated by different value^] of t, i.e. by 
different points p = 

Then we have the following result. 

Lemma 22. A is a finite set. 

Proof. It is clear that there are finitely many y 0 s satisfying (1), (2), (3) and 
(4). So let us see that there are also finitely many y 0 s satisfying (5). For this 
purpose, note that Od(C) has finitely many local singularities. Furthermore, 
since O c fiC) is simple by hypothesis and C is properly parametrized, then 
there are finitely many points of Od(C) generated by different values of the 
paramemeter t. Therefore, there are also finitely many yos satisfying (5). D: 


6 Recall that since Od(C) is irreducible by hypothesis, then its implicit equation consists 
of just one irreducible factor 

7 Notice that these last points correspond to self-intersections of Od{C) 
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Therefore, a generic y 0 does not satify any condition (l)-(5). This is 
crucial in the next result. 


Proposition 23. Let H(x,y ) = F(x,y) ■ G(x,y). Then F(x,y ) is a square- 
free polynomial. 


Proof. Since F(x,y ) = ( f(x,y )) r , we must prove that r = 1. In order to do 
this, observe first that f(x, y) cannot be of the form a-y, with «6l. Indeed, 
in that case C consists of a pair of parallel lines, which does not correspond 
to a rational curve. Therefore in order to prove the assertion, it suffices to 
prove that for a generic yo, F(x, yo) is square-free. For this purpose let y = yo 
be generic; hence, y 0 does not satisfy any condition (l)-(5). In particular, 
the leading coefficients of P(x,y,t ) and Q(x,y,t ) do not identically vanish 
for y = y 0 , so Res t (P(x,y,t),Q(x,y,t)) specializes properly, i.e. H(x,y 0 ) = 
Res t(P(x, yo, t), Q{%, 2/o, t)) (see Theorem [7] in Section 2.2). Since yo does not 
satisfy condition (1), the line y = y 0 does not intersect H(x,y ) = 0 in any 
point both belonging to Od(C) and to the curve G(x,y) = 0. Moreover y 0 
does not satisfy condition (4) either, and therefore we have that 


H(x, y 0 ) = (x- xx ) r ■■■{x~ x n ) r ■ G(x, y 0 ), 


where G(xi,y 0 ) ^ 0 for i = 1In other words, the intersection of 
the line y = yo with the curve F(x,y ) = 0 consists of the points (ojj, yo), 
i = 1 ,n, and all these points have the same multiplicity of intersection 
with y = yo, namely r. Now assume that r > 1. From Proposition 5 in [9], 
we have the following possibilities: 

(i) For any i = 1 ,...,n there is yi y yo and fo,^i with P(xj,yo,to) = 
Q(xi,y 0 ,t 0 ) = 0 and P(xi,y 0 ,ti) = Q(xi,y 0 ,ti) = 0. But this implies 
that ( Xi,y 0 ) G Od{C) is simultaneously generated by t 0 and t\, which 
cannot happen because y 0 does not satisfy condition (5). 


(ii) There is some i = 1,... ,n such that P(xi,y 0 ,t ) and Q(xi,y 0 ,t ) share 
a root to of multiplicity at least 2. However, in that case Pt vanishes 
at the point (ay,yo,to), which by Lemma [8] implies that (xi,yo) is a 
local singularity of Od(C). But this cannot happen because yo does not 
satisfy condition (5). 


(iii) The line x = x* is a common vertical asymptote of the curves (defined 
on the xt-plane) P(x, yo, t) = 0, Q(x, y 0 , t ) = 0. But this cannot happen 
because yo does not satisfy condition (3), and hence yo ^ y±oo- 
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So we conclude that r — 1. 


□ 
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